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Abstract 

fNj , We study the holographic QCD in the hadronic medium by using the soft wall model. We 

discuss the Hawking-Page transition between Reissner- Nordstrom AdS black hole and thermal 
charged AdS of which the geometries correspond to dcconfincmcnt and confinement phases 
respectively. We also present the numerical result of the vector and axial vector meson spectra 
r depending on the quark density. 
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1 Introduction 

Holographic QCD offers new approaches to understand the strongly interacting regime of gauge 
theories based on AdS/CFT duality pQ. The properties of QCD including confinement and chiral 
symmetry breaking have been discussed in asymptotically AdS geometries [21 HI [5j El [5] . 

The gravity dual of QCD can be derived in both top-down and bottom-up approaches. There 
are two methods to construct AdS/QCD background in the bottom-up approach. One is intro- 
ducing an infrared (IR) cutoff [H EJ [TOj [Til 021 El E], and the other is scaling the action or 
modifying the background metric by introducing a scalar field or a warp factor so that a smooth 
IR truncation is induced [151 fl6l fT7] . The light-front holography [18] is constructed by using the 
both methods. 

The phase transition between the confining and deconfining phases corresponds to the Hawking- 
Page transition [19] of gravity between thermal AdS (tAdS) at low temperature and the Schwarzschild 
AdS black hole at high temperature in pure gauge theory [2]. The phase transition of holographic 
QCD in the hard wall [UJ and soft wall model [15] was discussed in [20] . 

The holographic QCD can be improved by including the chemical potential [2JJ [22], E21 EH 
[251 [261 E3 [281 EH1 EQl EE]. The charge of Reissner-Nordstrom AdS black hole (RNAdS BH) 
can be interpreted as quark charge when branes fill the bulk [27]. The boundary value of the 
time-component gauge field is the chemical potential whose dual operator is the quark number 
density. The gravity dual to deconfining phase of QCD is the RNAdS BH. The thermal charged 
AdS (tcAdS) space, which is the zero mass limit of the RNAdS BH, was proposed as the dual 
geometry corresponding to the confining phase of QCD [28]. The Hawking-Page transition between 
confinement and deconfinement phases was studied in the presence of the chemical potential or 
the quark density in the hard wall model [281 EH]. It has been also observed that meson mass 
increases as the quark density increases in the hard wall model [30] . 

In this report, we discuss the holographic QCD in the hadronic (or quark) medium by using 
the soft wall model [15| where a non-dynamical scalar field is introduced to give a smooth IR 
truncation. The scalar field does not affect the dynamics of the gravity as it scales the full action. 
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In section 2, we investigate the Hawking-Page transition between RNAdS BH and tcAdS. In 
section 3, we study the vector and axial vector meson spectra depending on the quark density. In 
section 4, we summarize our results and discuss issues for future work. 



2 Hawking-Page transition in the soft wall model 

We consider the holographic QCD in the hadronic (or quark) medium using the soft wall model. 
According to the AdS/CFT correspondence, the quark number operator in the gauge theory side 
corresponds to the time-component of the bulk gauge field. So the dual geometry of the hadronic 
medium should be a five-dimensional spacetime containing some electric charges. There are well- 
known examples describing the asymptotic AdS space containing the electric charges. One is 
the Reissner-Nordstrom AdS black hole (RNAdS BH) corresponding to the deconfining phase 
described by quarks and gluons in the gauge theory side. The RNAdS BH in Euclidean spacetime 
signature is 

R 2 ( 1 \ 

ds 2 = -5- ( f RN (z)dt 2 E + - — --dz 2 + 5 ij dx i dx j , (2.1) 
z V Jrn(z) J 



with 



f RN {z) = l-mz A + q 2 z\ (2.2) 



where i,j imply the three-dimensional spatial indices. The parameters R, m and q are the AdS 
radius, the black hole mass and charge, respectively. We follow the convention of [29]. The 
time-component of the corresponding bulk gauge field is 

A = i (2tt 2 // - Qz 2 ) , (2.3) 

where /i and Q are related to the chemical potential and quark number density in the dual gauge 
theory. To satisfy the Einstein and Maxwell equations the quark number density Q should be 
related to the black hole charge q 

l3g 2 R 2 

In the confining phase the vacuum of the dual gauge theory is filled up with hadronic matter 
as quarks can not exist alone due to the confinement. The dual geometry of it was found in [28J, 
which is another solution satisfying Einstein and Maxwell equations. It is called a thermal charged 
AdS (tcAdS) space [28], whose metric in Euclidean signature is 

ds 2 = % (ftc(z)dt 2 E + j^dz 2 + Sijdx'dx?) , (2.5) 

Z \ Jtc{Z) J 

with 

ftc(z) = 1 + q' 2 z 6 . (2.6) 



Q=\l J h5-9- ( 2 - 4 ) 
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The time-component of the bulk gauge field is also in the form of (|2.3p . As will be shown, although 
this tcAdS solution is singular at z = oo it does not give any problem to calculate the physical 
quantities due to the potential barrier caused by the non-dynamical scalar field in the soft wall 
model. 

The Hawking-Page transition in the hadronic medium by using the hard wall model has been 
investigated [28|, [29] . The Hawking-Page transition in the soft wall model is also studied for the 
case of q = |15[ I20j. We investigate the Hawking-Page transition between the RNAdS BH and 
tcAdS backgrounds in the soft wall model, which describes the deconfinement phase transition of 
QCD in the hadronic medium. On these backgrounds, the Euclidean gravity action of the soft 
wall model is given by 



where the non-dynamical scalar field is $ = cz 2 . According to the AdS/CFT correspondence, the 
on-shell gravity action corresponds to the thermodynamic energy of the dual gauge theory, so we 
evaluate the on-shell gravity action on (|2.ip and (|2.5p with the gauge field (j2.3|) . 

For the deconfining phase, which is dual to the RNAdS BH, we impose the Dirichlet boundary 
condition Aq = i2-K 2 (i at the boundary z = to get the corresponding on-shell gravity action 



where V3, e and (3rn are the volume of the three-dimensional space, the UV cut-off and the inverse 
Hawking temperature /3rn = 1/Trn. 

As shown in [28} 129] . this on-shell gravity action is proportional to the grand potential of the 
dual gauge theory, which is a function of [1. The black hole charge should therefore be a function of 
the chemical potential \x. To determine this relation, we impose the regularity condition, Aq = 0, 
at the black hole horizon. Then, the black hole charge or the quark number density can be 
rewritten as 




(2.7) 




(2.8) 




(2.9) 



where z+ means the black hole horizon. The on-shell gravity action becomes 




(2.10) 



where 




(2.11) 
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with the exponential integral function Ei(z) 

EI(Z) : 



oo e -t 

dt . 

t 



(2.12) 



In ()2.10p . F(e) diverges in the limit as e — > 0, whereas the other terms remain finite. We renor- 
malize the on-shell gravity action (|2.10j) to have a well-defined grand potential of the boundary 
theory. We subtract the action of thermal AdS (tAdS or Euclidean AdS space) corresponding to 
the reference geometry, whose on-shell action is 



SAdS 



R 3 V H f 130 



R 3 V 3 



dt. 



00 4e~ 
dz = — 



(2.13) 



where /3q is the time periodicity in tAdS. F(oo) = has been used. Matching the circumferences 
between the RNAdS BH and tAdS at the UV cut-off e, (3q can be rewritten as 



A) 



v 7 Jrn{( 



Tbn 



Trn 



1 



m 



e 4 )+0(e b 



Using the expansion form of F(e) 



F{e) = -- i +0 



the renormalized on-shell action Srn in the e — > limit becomes 



Srn 



Srn - SAdS 
R 3 Vs 1 
k 2 T RN 



m 



F(z + ) + -+H{z + )-^- 



(2.14) 



(2.15) 



(2.16) 



As (I2.2p satisfies /(^+) = 0, we get the black hole mass m as a function of the black hole horizon 



m 



1 1 2 2 
— +qz + 



(2.17) 



Furthermore, the black hole horizon z+ can be expressed in terms of the Hawking temperature 
and chemical potential as follows 



3g 2 R 2 
8tt 4 kV 



l ^T 2 RN+ l ^£-,T RN \. 



3g 2 R 2 



(2.18) 
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Using the relations (|2.9p and (|2.17p . the renormalized on-shell action (|2. 16[) in terms of z+ and fi 
becomes 



e R S V 3 l 



T RN 



1 4tt 4 acV 87r 4 K V 2 



+ 24 + 3<? 2 # 2 4 + 3 5 2 i? 2 c4 ( e " + 



(2.19) 



Now, we calculate the on-shell action of tcAdS, whose dual theory lies in the confining phase. 
Substituting the solutions of (|2.3p and (|2.5p into (|2.7p . the on-shell action is obtained as follows 



S tc = ^ [' dt E J™ dz e- 2 (1 - 2</ 2 z) , 

where f3' is the periodicity in the Euclidean time direction. The integral range over z runs from 
the UV cut-off e to 00 because there is no black hole horizon for tcAdS. After performing the 
integration, the on-shell action for tcAdS becomes 



R 3 ^ Ql ( 1 



Stc = ^0 ( -F(e) - -e-^q' 2 ) , (2.20) 



c 



where F(oo) = and H(oo) = have been used. We renormalize the on-shell action due to the 
divergence of -F(e), as it is done for RNAdS BH. The circumferences of tcAdS and AdS at the UV 
cut-off should be the same. The relation between the Euclidean time periodicities is therefore 

Po = VMej f3' = P' + 0(e G ). (2.21) 

By subtracting the tAdS on-shell action (|2.13p . the renormalized on-shell action of tcAdS, St c in 
the limit as e — > becomes 



R 3 V ? 



i'2 



Stc = S tc - SAds = o~ $ — • (2.22) 



c 



The grand potential of the dual gauge theory is 

Vtc = y; (2.23) 

which is a function of the chemical potential. So q' should be represented as a function of [i. To 
determine q' as a function of fi, we first assume a simple relation between q' and \x 

q = an, (2.24) 

where a is a constant. Using this, the total quark number is obtained by the thermodynamic 
relation 



dn tc _ R 3 V 3 2aq' 

djJL K 2 



N = _<^l = il^±t^L. (2.25) 
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In [28] . it was shown that imposing the Neumann boundary condition instead of the Dirichlet one 
at the UV cut-off corresponds to the Legendre transformation from the grand potential to the 
free energy and the boundary action Sb is the same as [iN/3 1 . So, through the calculation of the 
boundary action we can determine the unknown parameter a. The boundary action at the UV 
cut-off is 



~ I ZiV'gW n M A^F 1 
9 JdM 



LN 

MN 9 



= ^-AapRyQ, (2.26) 
9 l 

where the unit normal vector is n = 

{o,0,0,0,-£v7tc(*)}- Comparing with the 

constant a is determined as 

a = 3c,r2 \W < 2 - 27) 

Substituting (|2.27p into (|2.22p . the renormalized action for tcAdS becomes 

B* = -**t*!t = J>!**V#f. (2.28) 
k c g z 

Usually the Hawking-Page transition occurs when the two on-shell actions for the RNAdS BH and 
tcAdS are the same. Before describing it, we first clarify (3' . Following [2\, f3' can be expressed 
in terms of the Hawking temperature Trjv by matching the time circumferences of these two 
backgrounds at the UV cut-off 



ft 



' fRN{e) 1 
/tc(e) T R n 

/ l-^ e 4 )+0(6 6 ). (2.29) 



Trn \ 2 

From this, the action difference between the on-shell actions for the RNAdS BH and tcAdS, which 
is proportional to the grand potential difference, becomes 



A5 = Srn — S, 



tc 



c 2 



^ + >-^>+^ + )-^-*' 2 > 



By expanding the F(e) as 



G 



the on-shell action difference becomes 
R 3 V 3 l 



AS 



V c2 + 



(cz 



1) + c z Ei(-czi) + 



1 



,2„2 



-cz± 9 



+ 



+ 



k 2 T RN \ z% K + y v +/ 24 ' 2 

where 2+, g and q' are functions of Trtv and ^ as shown in (|2. 18j) and below 



a? + si 

c c 



1 2 k 2 2n 2 fi _ 2N f 2tt 2 m 



We follow the convention of [27 



'3 K 2 „ 2 

2^ 2?r " C 



'3^V/ 2 
2ivf 27r ^ 



(2.31) 
(2.32) 



k 2 47r 2 i? 3 ' g 2 4^ 2 iT 



(2.33) 



for the last equalities. 

The Hawking-Page transition, which corresponds to the confinement/deconfinement phase 
transition occurs at AS = 0. We plot the confinement/deconfinement phase diagram for Nf/N c = 
0, 1/3, 2/3, 3/3 in Figureffl For N f /N c = 1 we find that n = 45.6MeV at T c = 0. This value is 
close to our expectation for the hadronic chemical potential. 

T [GeV] 

0.25 




"-- N f IN e 


= 


N f /N c 


= 1/3 


N f IN c 


= 2/3 


N f IN c 


= 3/3 



0.06 



0.08 



0.10 



H [GeV] 



Figure 1: The deconfinement temperature depending on the chemical potential. 
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3 Meson spectra in the soft wall model 



In this section, we investigate meson masses depending on the quark number density in the soft 
wall model. To describe the light mesons, we concentrate on tcAdS because it corresponds to the 
confining phase. Since we will not study the thermodynamic properties of the meson it is more 
convenient to use the Lorentzian tcAdS metric 

ds 2 = % (-ftc(z)dt 2 + j^dz 2 + SijdatdJ) , (3.1) 

z V Jtc{ z ) J 

where f tc (z) is 

f tc (z) = l + q 2 z 6 , (3.2) 

as defined in Section [2j The background gauge field corresponding to the quark number density 
is 

A) = 2vr 2 /i - Qz 2 , (3.3) 

in the Lorentzian signature. 

In the soft wall model [15], the gravity action describing meson spectra of the dual gauge 
theory is given by 



d 5 xV Z G e" $ Tr 



(3.4) 



where m 2 = — and the real part of <f) corresponds to the chiral condensate of the dual gauge 
theory. Tr means the trace over the flavor symmetry group indices. The superscripts, (L) and 
(R), in the last term imply the left and right parts of SU(Nf)i x SU(Nf)ji flavor symmetry 
group with /mn = d M a^" fl) — d N a[ L / R> — i \a^ ,K) } (% ' R) ~\ , where a M implies the bulk fluctuation 
corresponding to vector or axial-vector mesons. Usually, there exists non-zero chiral condensate 
in the confining phase which breaks the chiral symmetry spontaneously, so we set (j) as 

V ( Z ) i A2ir a T 



2\/2 



le 1 ^ 1 , (3.5) 



where v is real. 7r a will be considered as fluctuations corresponding to pseudoscalar mesons. After 
turning off all bulk fluctuations a M and 7r a , the equation of motion describing the chiral condensate 
is given by 

° = dz e ~* 9zV ) + ^ e ~" V > (3 - 6) 

where the gravitational backreaction of v is ignored and v is considered as a background field 
describing the chiral condensate like the original one [11]. Near the boundary z = 0, the asymptotic 
form of v is 

v = m q z + az 3 + • • • , (3-7) 
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where m q and a are the quark mass and chiral condensate respectively. The existence of the 
analytic solution of (|3.6p is not guaranteed so we will use the numerical solution depending on two 
initial parameters to investigate the meson spectra. We investigate the case of Nf = 2 and N c = 3. 
We set the five-dimensional coupling constant in (|3,4p . g 2 = 12 ^ R , following the convention of 



3.1 Vector meson 



We transform the vector fields and a^ 1 to obtain vector and axial vector fields 



(aS?+aJ?) 



2 

- (a (L) 
2 



) 



The action describing the vector and axial vector mesons up to the quadratic order is 



AS 



dPxV^Ge 



v 2 (d M 7T-a M )(d M 7r-a M ) + ^[ 



AV) 2 , f(A)2 

J MN ~T J MN 



(3.8) 



(3.9) 



where = d M v N 



d N v M and f[ A) 



i mn — d M a N — d N a M . We study the model in the gauge where 
v z = and a z = 0. Since Lorentz symmetry is broken we consider the vector field fluctuation in 
the 2,-direction. The Fourier transform of the vector field can be defined as 



Vi(t,z) = J ^Le-^Viiz), 



(3.10) 



with a condition rriy = to 2 . The vector field v M in the action (|3.9p does not couple to the scalar 
field v. The equation of motion for the vector field v% under the transformation (|3.10p is 



= d z ( &e-*d g vA + ^e"*7 
V z Zl ) Z f tc 

We introduce a dimensionless variable as follows 

2, = \[cZ. 

All the other dimensionful parameters get scaled accordingly 

^2-6 ~_ 1 



f tc (~z) = 1 + ? ~z b 



.3/2 ! 



(3.11) 



(3.12) 



(3.13) 



By the field redefinition v% = e x ipi with X = z 2 — log^, the equation (|3.1ip reduces to the 
Schrodinger-like equation 



X" X' 2 m 2 



(3.14) 



9 



Notice that the effective potential V e ff = — ( g — V nas the positive infinite value at 

two boundaries z = and z = oo, so the wave function ipi should be zero at these two boundaries. 
This implies that the natural boundary conditions for ipi are 

^(0) = Vi(oo) = 0. (3.15) 

In the case of zero quark density (q = 0, f tc = 1), the effective potential becomes V e ff = z 2 + 
^p- — fh v and the equation (|3.14p is exactly solvable. The solution for the differential equation is 

4^ = ^e-^LU~z% (3.16) 

with associated Laguerre polynomials L^x) = dL ^ x ^ ■ This is the (n + l)-th excited solution of 
(|3.14p for the quantized mass values rhy = m l l = 4(n+l) and corresponds to the (n + l)-th excited 
vector meson. From this result, we can see the Regge behavior of the vector meson, ~ n + 1 , 
in the zero density case. 

There is no analytic solution for the case of non-zero quark density. We will investigate 
the meson spectra numerically. It is worth noting that the solution (|3.16p satisfies the set of 
the boundary conditions (|3. 15[) . It shows that these are consistent boundary conditions for the 
numerical analysis. 

In Figure [2l we plot vector meson masses for n + 1 = 1, 2, 3, 4 depending on the quark density. 
The meson mass increases as quark density increases. This is qualitatively consistent with the 
vector meson spectra in the hard wall model [30J. 



q(GeV 6 ) 


n = 


n = 1 


n = 2 


n = 3 


0.0 x (0.388) 3 


0.776 


1.09743 


1.34407 


1.55200 


0.1 x (0.388) 3 


0.82285 


1.34296 


1.88084 


2.42465 


0.2 x (0.388) 3 


0.87629 


1.53163 


2.21219 


2.89822 


0.3 x (0.388) 3 


0.92270 


1.67681 


2.45737 


3.24278 


0.4 x (0.388) 3 


0.96361 


1.79728 


2.65726 


3.52160 


0.5 x (0.388) 3 


1.00035 


1.90151 


2.82844 


3.75933 


0.6 x (0.388) 3 


1.03386 


1.99413 


2.97948 


3.96848 


0.7 x (0.388) 3 


1.06477 


2.07793 


3.11547 


4.15639 


0.8 x (0.388) 3 


1.09274 


2.15479 


3.23970 


4.32777 


0.9 x (0.388) 3 


1.12037 


2.22599 


3.35445 


4.48587 


1.0 x (0.388) 3 


1.14602 


2.29249 


3.46134 


4.63300 



Table 1: The first four excitations of the vector meson spectra in GeV. yfc = 0.388GeV. 
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mass 2 [GeV 2 ] 




3.2 Axial vector meson 

We study the axial vector meson spectra. The axial vector can be decomposed into a transverse 
component and a longitudinal component as 

% = a M + d M X- (3-17) 

The a M corresponds to the axial-vector meson. We choose the axial gauge, a z = 0. We consider 
the fluctuation fixing a = as the Lorentz boost symmetry is not manifest. As the axial vector 
couples to the scalar field v in the action (j3.9[) we solve the equation of motion for the scalar field 
as well. Under the dimension scale (I3.12P the equation (I3.6P becomes 

d~ z (^e-" 2 d,v) + ^e-' z \ = 0. (3.18) 



The parameters in the scalar field ()3.7p are scaled as 

a 

c 3/2 • 



3 Q u / \ 

m q z + az H , m q = <J = -^- (3.19) 



By the Fourier transformed axial vector 



/ ^e-^%(co,z), (3.20) 



the equation of motion is obtained as 

ftc _<J> n - \ 1 



dz I ^ e -*cU; - -U - * R 2 9 2 v 2 - ^-m* 1^ = 0. (3.21 



~ A, l I Q \ J £ 

Z J Z 6 \ fi 



lu- 



ll 



We set R = 1 in (|3.2ip . By the field redefinition cij = e Y ^ 2 d with Y = z 2 — log the equation 
(|3.2ip can be rewritten as 

The mass spectra from the numerical solution of (|3.18p and (|3.2ip for m q = 0.005044GeV and 
a = (0.2619GeV) 3 are plotted in Fig [31 The meson mass increases as quark density increases. 



q(GeV 6 ) 


n = 


n = 1 


n = 2 


n = 3 


0.0 x (0.388) 3 


1.20195 


2.04792 


2.83813 


3.59508 


0.1 x (0.388) 3 


1.20415 


2.05663 


2.85662 


3.62591 


0.2 x (0.388) 3 


1.21044 


2.08104 


2.90742 


3.70920 


0.3 x (0.388) 3 


1.21996 


2.11703 


2.98027 


3.82575 


0.4 x (0.388) 3 


1.23155 


2.15995 


3.06494 


3.95831 


0.5 x (0.388) 3 


1.24391 


2.20579 


3.15386 


4.09546 


0.6 x (0.388) 3 


1.25581 


2.25161 


3.24237 


4.23097 


0.7 x (0.388) 3 


1.26636 


2.29591 


3.32850 


4.36251 


0.8 x (0.388) 3 


1.27537 


2.33861 


3.41222 


4.49003 


0.9 x (0.388) 3 


1.28327 


2.38044 


3.49422 


4.61415 


1.0 x (0.388) 3 


1.29080 


2.42216 


3.57504 


4.73521 



Table 2: The first four excitations of the axial meson spectra in GeV for m q = 0.005044GeV and 
a = (0.2619GeV) 3 . ^/c = 0.388GeV. 



4 Discussion 

We have studied the holographic QCD in hadronic medium by using the soft wall model, where the 
confinement scale is induced by a non-dynamical scalar field. We have observed the Hawking-Page 
transition between Reissner-Nordstrom AdS black hole and thermal charged AdS and obtained the 
phase diagram, which is analogous to the diagrams of the hard wall model. We have studied the 
vector and axial vector meson spectra depending on the quark density. In both cases the meson 
mass increases as quark density increases. The Regge behavior is not observed in the presence of 
the quark density. 

We have discussed the model in the backgrounds of RNAdS BH and tcAdS whose geometries 
are not modified by the scalar field which induces the IR cutoff. It would be interesting to study 
relevant physical quantities in the geometry whose metric gets deformed by a scalar field. We will 
report those results elsewhere. 
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mass 2 [GeV 2 ] 

20 - 




4th excitation 



15 - 



3rd excitation 



10 - 




2nd excitation 



1 st excitation 



o ^ 

0.00 



0.01 



0.02 



0.03 



0.04 



0.05 



0.06 



<7[GeV 3 ] 



Figure 3: The first four excitations of the axial meson spectra for m q = 0.005044GeV and a = 
(0.2619GeV) 3 . 
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